Update on tests of the Cen A neutron-emission model of highest energy cosmic rays 
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We propose that neutron emission from Cen A dominates the cosmic ray sky at the high end of the 
spectrum. Neutrons that decay generate proton diffusion fronts, whereas those that survive decay 
produce an angular spike in the direction of the source. We use recent data reported by the Pierre 
Auger Collaboration to normalize the injection spectrum and estimate the required luminosity in 
cosmic rays. We find that such a luminosity, Lcr ~ 5 x 10 40 erg/s, is comfortably smaller than the 
bolometric luminosity of Cen A, Lboi ~ 10 43 erg/s. We compute the incoming current flux density 
as viewed by an observer on Earth, and we show that the anisotropy amplitude is in agreement with 
data at the lcr level. Regardless of the underlying source model, our results indicate that after a 
decade of data taking the Pierre Auger Observatory will be able to test our proposal. 



About a decade ago we put forward the idea of us- 
ing neutrons as markers of ultrahigh energy cosmic ray 
(UHECR) emission from Cen A pj . In this Brief Report 
we update our proposal to accommodate recent observa- 
tions. 

The HiRcs Collaboration has reported a suppression of 
the CR flux above £ S HR = 56 ± 5(stat) ± 9(syst) EeV 0. 
The spectral index 7 of the flux, J oc E r , steepens 
from 2.81 ± 0.03 to 5.1 ± 0.7. This suppression has been 
confirmed by the Pierre Auger Collaboration, measuring 
7 = 2.69±0.2(stat)±0.06(syst) and 7 = 4.2±0.4(stat)± 
0.06(syst) below and above E^ = 40 EeV, respectively 
(the systematic uncertainty in the energy determination 
is estimated as 22%) 0]. In addition, an intriguing CR 
excess was found in the direction towards Cen A, a pow- 
erful radiogalaxy at d = 3.4 Mpc Out of the 69 
Auger events with E\ > 55 EeV (collected over 6 yr but 
equivalent to 2.9 yr of the nominal exposure/yr of the 
full Auger), the overdensity with largest significance is 
13 CRs within 18° from Cen A, versus only 3.2 expected 
if the flux were isotropic. Furthermore, 2 of the events 
arrived within less than 3° of the radiogalaxy. Within 
errors, we set £7™ ~ E^ ~ E\, and take the spectrum 
to fall approximately as J oc £V 4 above the onset of the 
suppression. 

HiRes has presented evidence that the CR composi- 
tion remains protons up to the highest energies Q . Auger 
data on the depth of shower maximum A" max , its rms fluc- 
tuation cr(V max ), and muon rates at ground level favor a 
heavy composition like 56 Fe Q . A critical assumption of 
our hypothesis is that the primaries of the highest energy 
Auger data are protons (and neutrons) and not 56 Fe. We 
note that some arguments have appeared recently which 
lend support to the hypothesis of proton primaries. One 
is the reasoning Q that if particles responsible for the 
Cen A excess are heavy nuclei (say 56 Fe), then a similar 
(as yet unobserved Q) anisotropy with better statistics 
should be present among protons at E\ /26 ~ 2 EeV en- 
ergies. A second is that modeling reveals that X max and 
cr(X max ) in themselves may be poor estimators of pri- 
mary composition 0. In [l(| it is shown that the ratio 
(k = A/ A) of the measured shower attenuation length (A) 



to the interaction length of protons in the atmosphere (A) 
is likely poorly estimated in shower codes, and poorly as- 
signed in various CR experiments. The assumption of a 
single value 1.26 ± 0.03 |10| for k, and proton primaries, 
is shown to give excellent agreement with all CR data 
as of 2007. The rate and fluctuations of early shower 
development (related to depth of first interaction X\ , in- 
elasticity K, mean multiplicity (n), and depth of most 
inelastic interaction X n ) are the physics behind the devi- 
ation of k from unity. In pd| it is shown that the simple 
combination AT max — <r{X max ) is a superior estimator of 
composition than V max or cj(V max ) individually. The 
combination is much less sensitive to data features that 
can skew interpretation, such as fluctuations in X± and 
tails in the V max distribution. Furthermore, when this 
new estimator is applied to the HiRes and Auger data, 
both data sets find agreement with a primary spectrum 
of protons [Tlj |. Thus, we feel that the assumption herein 
that UHECR primaries hitting the Earth's atmosphere 
are protons and neutrons is viable. 

We propose that neutron emission from Cen A dom- 
inates the observed CR flux above the suppression [1[. 
(The acceleration process of the parent proton popula- 
tion is discussed in Appendix [SJ) Neutrons that decay 
generate the proton diffusion fronts. The Bohm diffu- 
sive regime is usually described by an energy-dependent 
diffusion coefficient, D = 0.1 (E EcV /B nG ) (Mpc 2 /Myr). 
Charged particles with E > E c = diApc -BnG EeV propa- 
gate along a straight lines, whereas particles with E < E c 
diffuse [12]. For an extragalactic magnetic field B n Q = 50 
(see Appendix [Bj) , protons even up to the highest ob- 
served energies undergo Bohm diffusion. As a result of 
this diffusion, an injection spectrum of dNo/dEdt oc E~ 3 
in the region of the source cutoff, results in a spectrum 
J oc E~ 4 for the protons at Earth. The rate at Earth for 
the surviving neutrons is 



dNn 

dt 



S 



And 2 



Ei 



- d /X(E) 

dEdt 



(1) 



where X(E) ~ (£'/10 20 eV) Mpc is the neutron de- 
cay length and S is the area of the surface detector 
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(= 3000 km 2 for Auger). For the energy interval be- 
tween Ei = 55 EeV and E 2 = 150 EeV, we calculate the 
normalization factor using the observation of 2 neutrons 
in 3 yr. We then use this normalization factor to cal- 
culate the luminosity of the source in the above energy 
interval. We find = 9 x 10 39 crg/s. Next, we 

assume continuity of the spectrum at E\ as it flattens 
at lower energy to E~ 2 . Taking the lower bound on the 
energy to be E = 1 EeV, we can then fix the luminos- 
ity for this interval and find Lq£' Ei ' > = 4 x 10 40 erg/s. 
Adding these, we find the (quasi) bolometric luminosity 
to be L^''^ 2 ' = 5 x 10 40 erg/s, which is about a factor of 
2 smaller than the observed luminosity in 7-rays, L 7 
10 41 erg/s QJ] in the interval 100 McV < E < 10 GeV. 

To further constrain the parameters of the model, 
we evaluate the energy- weighted approximately isotropic 
(diffuse) proton flux at 70 EeV. The value of this diffuse 
flux depends on the energy threshold E\ above which Cen 
A is the dominate source (55 EeV here), on the spectral 
index at the source (—3 here), on the nature of diffusion 
(Bohm here), and on the on-time T on of the source [l[. 
We obtain (E 4 J) = 1.6 x 10 57 eV 3 km -2 yr" 1 sr" 1 , in 
agreement with observations Q, for a source actively 
emitting UHECRs over T on > 70 Myr. (We note that 
the minimum T on is similar to the diffusion time Tjj (E) ~ 
d 2 /[4D(E)} ~ 20 Myr at E = 70 EeV and B nG = 50.) 
If we assume circular pixel sizes with 3° radii, the neu- 
trons will be collected in a pixel representing a solid angle 
AO ~ 8.6 x 10~ 3 sr. The proton rate coming from the 
direction of Cen A is 



dN p 
~~dT 



r E2 dE 
= S AQ (E J) —r = 0.08 events/yr 
Je x e 



(2) 



It is important to stress that the 3° window does not 
have an underlying theoretical motivation. Recall that 
this angular range resulted from a scan of parameters 
to maximize Auger's signal significance. Cen A covers 
an elliptical region spanning about 10° along the major 
axis. Therefore, some care is required to select the region 
of the sky which is most likely to maximize the signal-to- 
noise jl4l | . 

We now address the question of anisotropy. Duplicat- 
ing the analysis of |l[ for Bohm diffusion on a 50 nG B- 
field we determine the amplitude of a dipolc term aligned 
with Cen A, the so-called "asymmetry parameter" a. For 
E = 70 EeV we obtain a = 0.29, within la of the same 
anisotropy amplitude a = 0.25 ± 0.18 obtained from the 
arrival directions of the 69 observed events [l5| . 

One caveat is that we assumed that neutrons com- 
pletely dominate the ultrahigh energy Cen A emission 
spectrum; that is 



dN 
dEdt 



= (N£ + NP)E~ 3 , with iVj/W n <l. (3) 



This reduces the number of free parameters in the model. 
The actual proton-to-neutron fraction depends on the 
properties of the source, especially the ratio of photon- 
to-magnetic energy density. 



In summary, existing data is consistent with the hy- 
pothesis that Cen A dominates the CR sky at the high 
end of the spectrum [l6| . The Pierre Auger Observatory 
is in a gifted position to explore Cen A and will provide in 
the next 9 yr of operation sufficient statistics to test this 
hypothesis. We expect about 6 additional direct neutron 
events against an almost negligible background fl7j . 

The potential detection of neutrons at Auger can sub- 
sequently be validated by the larger aperture of space- 
based UHECR experiments. The JEM-EUSO mission 
is scheduled to launch in 2017, and remain operational 
aboard the International Space Station for 3 or 5 years. 
Including the 20% duty cycle and 50% Southern ver- 
sus Northern hemisphere exposure, one expects 60-100 
events within 18° of Cen A for the three- year JEM-EUSO 
mission, and therefore 9-15 direct neutrons. For the five- 
year mission, numbers are proportionately higher. Work 
on the implications of our Cen A model for the JEM- 
EUSO mission is in progress fl4| . 
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Appendix A: UHECR emission from Cen A 

Cen A is a complex radio-loud source identified at opti- 
cal frequencies with the galaxy NGC 5128 [HI]. Radio ob- 
servations at different wavelengths have revealed a rather 
complex morphology. It comprises a compact core, a jet 
(with subluminal proper motion /3j Gt <~ 0.5 [l9j]) also vis- 
ible at A"-ray frequencies, a weak counter-jet, two inner 
lobes, a kpc-scale middle lobe, and two giant outer lobes. 
The jet would be responsible for the formation of the 
northern inner and middle lobes when interacting with 
the interstellar and intergalactic media, respectively. 

In order to ascertain the capability of Cen A to ac- 
celerate UHECR protons one first applies the Hillas cri- 
terion [20l | for localizing the Fermi engine in space, i.e. 
that the Larmor radius be less than the size of the mag- 
netic region. From this condition we infer a maximum 
CR energy of 



E — B^q i?k pc EeV . 



(Al) 



Furthermore, since the magnetic field carries with it an 
energy density B 2 /8ir and the flow carries with it an 
energy flux > f3cB 2 /(8tt), (|A1[) also sets a lower limit on 
the rate 
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Lb > g /3jctci? kpc -u MG 



~B„ 



(A2) 
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at which the energy is carried by the out-flowing plasma, 
and which must be provided by the source [2l| . The min- 
imum total power of the jets inflating the gian t lobes of 
Cen A is estimated to be « 8 x 10 43 erg/s |22|]. This ar- 
gument provides a conservative upper limit for the mag- 
netic field in the jet with kpc-scale radius, B m g ^ 50, 
and through (fM) leads to E = 50 EeV H]. 

Of particular interest here, it was recently noted that 
shear acceleration [24| could help push proton energies up 
to and beyond 50 EeV [25[. The limb-brightening in the 
X-ray jet together with the longitudinal magnetic field 
polarization in the large scale jet might be indicative of 
internal jet stratification, i.e. a fast spine surrounded by 
slower moving layers. Energetic particles scattered across 
such a shear flow can sample the kinetic difference in the 
flow and will naturally experience an additional increase 
in energy. It is therefore of interest to explore the ex- 
treme case, in which protons that diffuse from the inner 
shock region into the outer shear layers charge-exchange 
to produce neutrons with E 2 ~ 150 EeV. (We let exper- 
iment be the arbiter of this proposed mechanism.) 

Appendix B: Extragalactic Magnetic Field 

Surprisingly little is actually known about the extra- 
galactic magnetic field strength. There are some mea- 
surements of diffuse radio emission from the bridge area 
between the Coma and Abell superclusters (26[, which 
under assumptions of cquipartition allows an estimate 
of 0(0.2 — 0.6) fiG for the magnetic field in this region. 
Fields of 0(/iG) are also indicated in a more extensive 
study of 16 low redshift clusters [27|- It is assumed that 
the observed B-fields result from the amplification of 
much weaker seed fields. However, the nature of the 
initial week seed fields is largely unknown. There are 
two broad classes of models for seed fields: cosmologi- 
cal models, in which the seed fields are produced in the 
early universe, and astrophysical models, in which the 
seed fields are generated by motions of the plasma in 
(proto)galaxies. Galactic winds are an example of the 
latter. If most galaxies lived through an active phase 
in their history, magnetized outflows from their jets and 
winds would efficiently pollute the extragalactic medium. 
The resulting i?-ficld is expected to be randomly oriented 
within cells of sizes below the mean separation between 
galaxies, As < 1 Mpc. 

Extremely weak unamplificd extragalactic magnetic 
fields have escaped detection up to now. Measurements 
of the Faraday rotation in the linearly polarized radio 
emission from distant quasars [28| and/or distortions of 
the spectrum and polarization properties in the cosmic 



microwave background (CMB) [2|||3(| imply upper limits 
on the extragalactic magnetic field strength as a function 
of the reversal scale. It is important to stress that Fara- 
day rotation measurements (RM) sample extragalactic 
magnetic fields of any origin (out to quasar distances), 
while the CMB analyses set limits only on primordial 
magnetic fields. The RM bounds depend significantly 
on assumptions about the electron density profile as a 
function of the redshift. When electron densities follow 
that of the Lyman-a forest, the average magnitude of the 
magnetic field receives an upper limit of B ~ 10 -9 G for 
reversals on the scale of the horizon, and B ~ 10 -8 G for 
reversal scales on the order of 1 Mpc [3l[. As a statis- 
tical average over the sky, an all pervading extragalactic 
magnetic field is constrained to be [HJ 

VO-02 J \0.72J \MpcJ ' V ' 

where H^h 2 ~ 0.02 is the baryon density and h ~ 0.72 
is the present day Hubble expansion rate in units of 
100 km/sec/Mpc. (This is a somewhat conservative 
bound because fih has contributions from neutrons as 
well as from protons, but only electrons in ionized gas 
are relevant to Faraday rotation.) 

Very recently it was pointed out that the study of 
the energy-energy-correlation (EEC) in a given sample 
of UHECRs can provide important clues on the all per- 
vading extragalactic magnetic field (33| . Since it is ex- 
pected that UHECRs are accelerated at discrete sources, 
the deflection in cosmic magnetic fields would result in 
an energy ordering in the distribution of arrival direc- 
tions. The Pierre Auger Collaboration has recently re- 
leased their analysis of EEC [H|. The measured EEC 
distribution is compatible with the expectation from 
isotropic arrival directions, i.e. no energy-ordered de- 
flections are observed near the most energetic UHECRs. 
Such an uncorrelated distribution can be caused either 
by a high source density for an isotropic source distri- 
bution or by large deflections of the UHECRs in cosmic 
magnetic fields. Since the nearby distribution of matter 
is anisotropic we conclude that the second option is more 
viable. One interpretation of Auger data on EEC uses 
the PARamctrizcd Simulation Engine for Cosmic rays 
(PARSEC), to obtain a lower bound of 10 nG at 95% 
CL on an all-pervading extragalactic magnetic field with 
coherence lengt h Ar = 1 Mpc, assuming a source de nsity 
< 10" 4 Mpc" 3 This agrees with the results of [Iff. 
Use of these numbers should be regarded as tentative 
while awaiting additional data and further independent 
analyses. 
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